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We study the quantum brachistochrone evolution for a system of two spins- 1 
describing by an anisotropic Heisenberg Hamiltonian without zx, zy interecting 
couplings in magnetic field directed along z-axis. This Hamiltonian realizes quan- 
tum evolution in two subspaces spanned by | tt)> I 44) an d | 44) > I 44) separately 
and allows to consider brachistochrone problem on each subspace separately. Using 
operator of evolution for this Hamiltonian we generate quantum gates, namely an 
entanler gate, SWAP gate, iSWAP gate. We also show that the time required for 
the generation of an entangler gate and iSWAP gate is minimal from all possible. 

PACS number: 03.65.Xp, 03.65.Aa, 03.65.Ca, 03.67.Lx 

1 Introduction 

The quantum brachistochrone problem is formulated similarly as a classical 
brachistochrone [H El El SI El E] : What is the optimal Hamiltonian, under 
a given set of constraints, such that the evolution from a given initial state 
l^j) to a given final one \ipf) is achieved in the shortest time? For the first 
time this problem was considered by Carlini et al. [TJ. Using the variational 
principle, they presented a general framework for finding the time of optimal 
evolution and the optimal Hamiltonian for a quantum system with a given set 
of initial and final states. The authors solved this problem for some specific 
examples of constraints. In [2] it was shown that analogous results as in [JJ 
can be obtained more directly using symmetry properties of the quantum 
state space. Their approach was based on the idea considered in [7], where 
an elementrary derivation was provided for the minimum time required to 
transform initial quantum state into another orthogonal state. In [8] the 
quantum brachistochrone problem for mixed states was considered. 

Bender et al. explored the brachistochrone problem for a PT-symmetric 
non-Hermitian Hamiltonian [3] and showed that for the non-Hermitian PT- 
symmetric Hamiltonians satisfying the same energy constraint the optimal 
time of evolution between the two fixed states can be made arbitrarily small. 
For a more detailed discussion on this subject see in jU |3j. The quantum 
brachistochrone problem for a spin-1 system in the magnetic field was solved 
in [9]. 
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In [TUj [TTj [T2"| [T5] it was established that the entanglement of quantum 
states is impotant in connection with optimization of quantum evolution. It 
was discovered by Giovannetti, Lloyd and Maccone [T0| [TT] that, in certain 
cases, entanglement enhances the speed of quantum evolution of composite 
systems. Conection between entanglement and time of evolution in the case 
of two-qubit and N-qubit systems was explored in [TJ]. The authors of [H] 
showed that, as the number of qubits increases, very little entanglement is 
needed to reach the quantum speed limit. Also, it was recently explored that 
entanglement is an essential resource to achieve the speed limit in the context 
of quantum brachistochrone problem [151 EH EE] • 

In [18] it was formulated a variational principle for finding the optimal 
time of realization of a target unitary operation, when the available Hamil- 
tonians are subject to certain constraints dictated by either experimental or 
theoretical conditions. This method was illustrated for the case of a two-spin 
system described by an anisotropic Heisenberg Hamiltonian with Ja interac- 
tion between spins (i = x, y, z) and magnetic field h" (a = 1, 2 for the first 
and second spin, respectively) which is directed along z-axis. Using operator 
of evolution for this Hamiltonian there was generated three examples of tar- 
get quantum gates, namely the swap of two qubits (Uswap), the quantum 
Fourier transformation and the entangler gate (Uent)- In this case, they 
showed that the optimal time required to generate these quantum gate are 
not the shortest possible time. 

It was shown that the two-qubit Hamiltonian with the J xx and J yy inter- 
action combines CNOT gate with the iSWAP operation [19]. The quantum 
CNOT gate is the fundamental two-qubit gate for quantum computation ]20] . 
This gate plays a central role in networks for quantum error correction [21]. 

In this paper we consider the quantum brachistochrone problem for the 
system of two-qubit represented by Heisenberg Hamiltonian with J u (i = 
x,y,z), Jjk (J 7^ k = x,y) interaction between spins and magnetic field h% 
(a = 1, 2 for the first and second spin, respectively) which is directed along z- 
axis. This Hamiltonian realizes quantum evolution in two subspases spanned 
by | tt)j I 11) an d I T4)) I It)- We solve this problem for each subspase 
separately and using operator of evolution for this Hamiltonian we construct 
entangler, SWAP and iSWAP gates. Contrary to [18], our Hamiltonian con- 
tain J xy and J yx couplings which allow to generate the entangler gate during 
the shortest possible time. Also, the generation time of the iSWAP gate is 
the minimal possible time which corresponds to the brachistochrone one. 

Thus, this problem might have an important application in quantum 
computing, quantum teleportation and quantum cryptography, because, as 
we are going to show, it allows to reach maximally entangled states during 
the shortest possible time. 
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The paper is organized as follows. In Section [2] we introduce and dis- 
cuss the Hamiltonian of a system of two spins- 1. In Section [3] we solve the 
quantum brachistochrone problem on the subspace spanned by | ft)> I Tl) 
(subsection 13. ip and on the subspace spanned by | tl) , I II) (subsection 13. 2p 
separately and obtain Hamiltonians which provide optimal evolution on each 
subspaces. We obtain optimal conditions and time required to generate an 
entangler, SWAP and iSWAP gates in Section HI Finally, the summary and 
discussion are given in Section HI 



2 Hamiltonian 

Our aim is to explore quantum brachistochrone evolution for a system of 
two-qubit. What we consider is a physical system of a two-qubit represented 
by two-spin interacting via anisotropic couplings J a (i = x,y,z), Jjk (j ^ 
k = x,y) and magnetic field h° (a = 1,2 for the first and second spin, 
respectively) which is directed along z-axis. In other words, we choose the 
following Hamiltonian: 



H 



E 

i,j=x,y,z 



T ^!^ 2 



E 

i 



a 



where J xz = J zx = J yz = J zy = 0, a\ = Oi <8> J, of = I <g) a i} Oi are the 
Pauli matrices. Note that this Hamiltonian does not contain items with 
(T x a 2 z , o\a\ and CyO"^, o^cXy, therefore, it realizes quantum evolution on two 
subspaces spanned by | ff), | H) and | fl); I It) an d does not mix these 
subspaces. This quality allows to rewrite Hamiltonian ([T|) as H = Hi + Hjj, 
where in the basis labeled as | tt)> I I It); I 11); the Hamiltonians 
read: 



Hi 



H 



n 



( h + + J zz 



V J Re + 
/ 

hr - J zz 

J Re ~~ i"Jlm 

V 



J Re i"Jtm \ 




-h + + J zz J 
0\ 

Re + iJlm 


/ 








-hr - J ; 
o 



(2) 



where we have introduced h = 
Hi and H u commute ([Hi, H n 



h\±h 2 z , J fa 



Jxx^ Jyy and J i m 



Jxy i Jyx • 



0). Accordingly, Hi realizes the evolution 



3 



of a system on the | ft) > I 44) subspace and H n realizes the evolution on the 
I t4)j I 4t) subspace. 

Hamiltonian ([I]) has four eigenvalues: Ef = J zz + ui, Ej = J zz — ooi, 
Ef = —J zz + Uu, Ejj = — J zz — un with the corresponding eigenvectors: 



1^7) = 

|V&> = 

Nfr) = 

where we introduce uj[ = \l J 



2uj(uji — h + ) 
1 

2cu I (u I + h+) 
1 

2cuii(ojii — h~ 
1 

2uj n {tjj n + h' 



(j R e - iJL) I tt) + fa - h + ) I 44)] , 
{j R e ~ iJfm) I tt) " (W, + h + ) | 44)] , 

(4e + iJF m ) 1 14) + (wj/ - hr) | 4t)] , 

{jRe + iJTm) I t4) " ("II + ^) I 4t)] , (3) 



'Re 



2 + Jfm + an d U H 



J Re ^Tm ~\~ h 



Hamiltonians Hi and Hjj have common set of eigenvectors ([3]) • Hamilto- 
nian Hi has two eigenvalues Ef, Ej with corresponding eigenvectors \ifrf), 
and one two-fold degenerate eigenvalue with \ipf) and \ipJi) eigen- 
vectors. Similar situation is in the case for Hn. It has two eigenvalues Ef T , 
Ejj with eigenvectors \ipfi), l^n), respectively and one two-fold degenerate 
eigenvalue with two eigenvectors \ipf) and 



3 Quantum brachistochrone 

The brachistochrone problem is as follows: What is the optimal Hamiltonian, 
with the finite energy condition, such that the evolution from a given initial 
state to a given final one \ipf) is achieved in the shortest possible time. 
The finite energy condition can be written in the form [IB] : 

Tr H 2 - 2cu 2 = 0, (4) 

where a; is a constant. The finite energy condition (J4]) for Hamiltonian ([1]) 
is: 

ui 2 + u n 2 + 2J zz 2 = u 2 . (5) 

As we mentioned earlier, Hamiltonian ([T]) realizes evolution in two sub- 
space separately and does not mix these subspaces due to which we cannot 
observe evolution between states from different subspaces. Consequently, we 
can consider quantum brachistochrone problem on each subspace separately. 
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Let us consider the operator of evolution with Hamiltonian ([I]): 

U{t) = e lHt = e iHlt e lHllt = Ufa, (6) 

where we use that Hj and Hjj commute, and Uj(t) = e tHlt , Ujj(t) = e lHnt . 
We take % = 1. Uj realizes the transformation = Uj\tpi) on the subspace 
spanned by | tt)) I II) an d acts as a unit operator on another subspace 
(subspace spanned by | tl)> | it))- Contrary to Uj, Un realizes similar 
transformation \ifif) = Uu\ipi) on the subspace spanned by | ||), I It) an d 
acts as a unit operator on the subspace spanned by | tt)j | Li). Let us 
consider quantum brachistochrone problem on each subspace in detail. 



3.1 Quantum evolution on the subspace spanned by 

in), ui) 

We consider the quantum brachistochrone problem of two spins- 1 represented 
by Hamiltonian ([1]) on the | tt)> I II) subspace. As we noticed earlier Hamil- 
tonian H ([[]) acts on any state = a\ ft) II) (where the normalization 
condition is the following: |a| 2 + |6| 2 = 1) which belongs to subspace spanned 
by | tt), I II) following: 

H\iP)=Hj\7P). 

Therefore, we consider the quantum brachistochrone problem on this sub- 
space using Hamiltonian Hi from ([2]). Let us introduce the following opera- 
tors: 
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(7) 



These operators satisfy properties of Pauli and unit matrices for | tt) ? I II) 
states when we denote their as follows: | tt) = I t)> I II) = I !)• 

With the help of introduced operators <^ the Hamiltonian Hi from (j5J) 
can be written in the form: 

Hi = ct 1 ■ h 1 + J ZZ I J , (8) 
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where h 7 = (J Re , Jj~ m , h + ). This Hamiltonian is similar to the Hamiltonian 
of one spin- 1 in the magnetic field. In the case of one spin-| vector, h 7 is 
vector of the magnetic field. Brachistochrone problem for the spin-| in the 
magnetic field was considered in the paper [6]. We cannot use the result from 
this paper because we assume another finite energy condition (J4]) comparring 
to paper [B|. In [H] author fixed the largest and the smalest eigenvalues of 
the Hamiltonian. 

Now, using ([8]) we represent the evolution operator Ui = e lHlt as follows: 



I-I 1 



where 
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Here we use that (& 1 ■ h 1 
n = 1, 2, 3, ... In a matrix form 



vJi""! 1 and (cr 1 ■ h 1 



2n 



2n+l 



u)i 2n a l ■ h 1 where 



( (cos (ujt) + jj- sin (to it) h + )e 



U, 



\ 






^sin (uiit) (Jr 6 + iJ/ m )e 
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u„t o 
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^sin (ujt) (Jfe — iJ't 



Im) 



iJzzt 






cos to it — ^ sin (to it) h + )e 



Let us put the initial state as = \ tt) and the final one as \tpf) = 
a\ tt) + l-D- Then using the matrix representation for the operator of 
evolution Ui ffTUl) . the relation \ipf) = Ui\ipi) takes the form: 



I a \ 
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V 4 Sin (uit)( J Re + iJ Im) J 



From the fourth component of ffTTl) we obtain the necessary condition that 
the initial state reaches the final one: 



LOl 



sin (u T t) V J Re +Jj r , 



\b\. 



A2) 



From equation (1121) we obtain the time required to transform the initial state 
| ft) i nt o the final one a\ tt) + b\ ||): t 



arcsm 



1 + 



Re 1 im 
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Now, using the finite energy condition (jSJ) we optimize the time of evolution 
if we put the following condition: 



oo n = J zz = h + = 0. (13) 

The optimal time is thus: 

r = — arcsin |6|. (14) 

OJ 

If the final state is \ipf) = \ 44) then we obtain the following passage time 
r = 2~ (The shortest time of evolution between the two fixed orthogonal 
states is called the passage time [7]). 

Condition (j!2p does not account the phase of b. Let us find the condition 
that allows to reach the phase of components of the final state during the 
optimal time (|14p . For this we insert the optimal time of evolution f lT4|) 
and conditions ( 1T3"|) in equation ffTTj) . From the fourth component of ffTTj) 
we obtain the following equation: ^jb\ [JRe+iJfr^J = b, which allows to find 
another conditions for optimal evolution: J^ e = Jf m = — w^, where 
b = 9% + i^sb. Hence, all the necessary conditions for the optimal evolution 
in the subspace spanned by | ft), | LL) read: 

oj ^sb oj $i.b 

Jxx — ~Jyy = 2"T^T' '^ xy = ^ yx = ~ ~2 T^T' ^ zz = ^' h z = h z = 0. (15) 



Initial state can reach maximally entangled states (| tt) + e% ^\ 44) 
[0, 2ir]) during the minimal time r = 

If we put optimal conditions (ITS"]) on the Hamiltonian H (JTJ we obtain 
Hamiltonian which provides optimal evolution on the subspace spanned by 
I tt)) I 44)- 111 the matrix form it reads: 



H 



( zcj^ \ 





V-^R J 



(16) 



The Hilbert space of the | tt); I 44) subspace is two-dimentional and is 
represented by the Bloch sphere (the Bloch sphere is the state space of a 
two- level system with the unit radius). Any pure state can be identified as a 
point on this sphere. The shortest distance between |^) and \if)f) is a large 
circle arc on the sphere (geodesic path). The optimal way of transporting 
into \i/jf) is therefore to rotate the sphere around the axis ortogonal 
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to the large circle. The axis of rotation passes along two quantum states 
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s/2 
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41} which are eigenvectors of optimal 
Hamiltonian H ( I16p that correspond to the largest u) and the smallest —u 
eigenvalues. 



3.2 Quantum evolution on the subspace spanned by 

in>, ut) 

Similar result we obtain for this case. We consider the quantum brachistrochrone 
problem using Hamiltonian Hjj from ([2]), because Hamiltonian H ([I]) acts on 
any state \vf)) = a\ fl) + b\ It) (where the normalization condition is as 
follows: \a\ 2 + \b\ 2 = 1) which belongs to subspace spanned by | f^), | |t) 
following: 

H\1,} = H II \iJ>). 

Let us introduse the following analogues Pauli and unit matrices for the basis 
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Using (fT7|) and making the same steps as in the previous case we rewrite 
Hamiltonian Hjj from ([2]) as: 
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h /7 - J ZZ I 
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the components of vector h 11 read h 11 = (J^ e , —Ji m , h ). The operator of 
evolution for this case Uu = e lHnt we can represent as: 
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Here we use that i^a 11 ■ h 7/ ) = u n 2n I n and (a 11 ■ h 
where n = 1, 2, 3, ... In the matrix representation Uu reads as: 



, , In— II ull 
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t4) and the final one as 



ti 



Let us put the initial state as 
) + b\ |t) • Similarly as in previous case using matrix representation for Uu 
f[2"Uj) . we write the relation l^j) = Uji\ipi) in the form: 
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In this case conditions for optimal evolution read: 
2 W 



•7, 



7ry 



J, 



2 W 



J« = 0, ^ = ^ = 0, (22) 



and the optimal time is (|T4|) . Also, in this case, the initial state | tl) can 
reach maximally entangled states (J tl) + e **| It)) (0 £ [0,27r]) during 



the minimal time t = 



If we put conditions (1221) on the Hamiltonian H ([T]) we obtain Hamiltonian 
which provides optimal evolution on the subspace spanned by | tl); I It)- m 
the matrix form this Hamiltonian reads: 
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(23) 



The Hilbert space of the | tl) ; I It) subspace is two-dimentional and 
it is represented by the Bloch sphere. In this case the axis of rotation 



passes along two quantum states ^ | tl) — £]fr| It) 
which are eigenvectors of the optimal Hamiltonian H 
the largest u and the smallest —u eigenvalues. 
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4 Realization of quantum gates by two inter- 
acting spins 

Operator of evolution can be related with some quantum gates. We now 
demonstrate this explicitly by a few examples. Let us consider the entangler 
gate: 



where (f) € [0, 7r]. If this gate acts on the state | "ft), it will produce the 
^-dependent entangled state cos<^| ft) — sin0| ||). The comparison of §6§ 
with fl24l) using ffTQ]) and fl20|) leads to the following set of parameters: J xx = 
Jyy = Jzz = 0, J xy = J yx = ~, h\ = h 2 z = and t = ^. If we compare 
t = - with optimal time f[T4"j) we can see that this time is the optimal time of 
evolution along the brachistochrone on the subspace spanned by | "tt)> I ID 
and the connection between parameter (p and b is as follows <p = arcsin \ b\ (4> 
is called the Wooters distance). 

Entangled gate allow to reach Bell states from the nonentangled states. 
When (ft = |, this allows reaching the maximally entangled Bell state |$ + ) = 
75 (I ft) + I ID) fr° m tne initial state = | ID an d Bel1 state = 
"75 (I ^) — I ^D) from the initial state | tt)- Time of evolution between 
these pairs of states is the shortest possible time t = j-. Initial state | ft 
) can reach ortogonal final states | \D) during the time t = which is 
the passage time. In [T5] the entangler gate is generated from operator of 
evolution which is obtained from Heisenberg Hamiltonian with Ja interaction 
between spins (i = x,y,z) and magnetic field (a = 1,2 for the first 
and second spin, respectively) which is directed along z-axis. In paper [18] 
the generation time of the entangler gate is not the shortest possible time, 
because this Hamiltonian does not contain items with J xy and J yx interaction 
which guarantee the shortest time of generation Uent- 

Let us construct a similar gate which produces the entangled states on the 
subspace spanned by | t|), | ]/[)■ This gate is constructed with the entangler 



/ cos sin cj) \ 



Uent 



10 
10 



(24) 
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V o o o 








1 / 



(25) 
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where U NOT = U^ot ® I- If gate (j5SJ) acts on the initial state | ||) we get 
entangled state cos0| tl) + sin0| |t). Similarly to the previous example 
we compare unitary operator ([6]) with quantum gate (|25|) and obtain the 



following set of parameters: J xx = J yy = J zz = 0, X 



h 



X y — Uy X — 2 , 

ENT is the optimal time of 



h 2 z = and t = -. Time of generation of U- 
evolution along the brachistochrone (fl4|) (0 = arcsin|&|) in the subspace 
spanned by | ||), I It)- 

Let us insert = | in U-^j^ (I25|) and act by Uj^p on the initial 
nonentangled state | ti) or I it)- We reach the maximally entangled Bell 
state |^+) = i (| t|) + | It)) or the state that discribe EPR pair |^ _ ) = 
75 (I tt) — I It))) respectively. In this case, the time of evolution between 
these states is also the shortest possible time t — Evolution between two 
ortogonal states | tl) an d | It) is realized within the passage time t — 

The next important gate which we consider is the SWAP gate: 



SWAP 



^ 1 ^ 

10 

10 

V o o o i J 



(26) 



which exchanges the states of two qubits. Similar steps as we did in case for 
entangler gate, we obtain the following set of parameters: J xx = J yy = J zz = 



and t = As we can see, 



J ( J — ±JiL) J = J — n h 1 = h 2 

such as in |18j . time of generation of SWAP gate is not the shortest time, 
but the evolution is along a geodesic path. This is because the generation 
speed v = % is less then the maximal value of speed for this Hamiltonian 

"Umax 2CJ. 

As in the previous example we can generate the gate which exchanges the 
| tt) state into the | ||) and vice versa: 
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NOT 



( 1 \ 

10 
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V i o o o J 



(27) 



when we put the following conditions: J x 
Jt.„ = J,™ = 0, hi = hi = 0, and t 



=C v / 6 J 



t is also 



'.).)■ Jyy Jzz J \ J 

>x V - "yx - ,« z - ,« z - w, » — ^ on unitary operator ©. 
not the shortest time, because the generation speed v = % is less than the 
maximal speed v max = 2cu, but evolution is also along a geodesic path in this 
subspace. 
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As a last example, we consider the unitary operator UiswAP'- 



U, 



iSWAP 



^ 1 ^ 

i 

i 

V 1 J 



(2f 



It can be obtained from the time unitary operator ([ 
conditions. J ra Jyy ~2> ^zz Jxy Jyx 0, h z 
Also, we can obtain the following gate: 



if we put the following 

hi = and t = 

z 2oj 



u- 



iSWAP 



NOT u iSWAPU NOT 



ul 



( i \ 

10 

10 

V i J 



(29) 



which is constructed using Uiswap and U^ot- This gate can be obtained 



if we put the following conditions: J xx = —J yy = tt, J zz = J xy = J yx = 0, 

hi - Hi 



and t = on unitary operator Q. 
The generate times Uiswap and U iSWAP are the shortest possible times 
of evolution between the ortogonal states on the subspaces spanned by | ti) , 
| and by \ft), \ 41), respectively 



5 Summary and Discussion 

We have studied the quantum brachistochrone problem for the system of 
two-qubits represented by two spins interacting via anisotropic couplings Ju 
(i = x,y,z), Jjk (j 7^ k = x,y) and magnetic fields h® (a = 1,2 for the 
first and second spin, respectively) which is directed along z-axis (pQ). This 
Hamiltonian realizes quantum evolution in two subspaces spanned by | f[), 
| If) and | f4), I It) and does not mix these subspace because it does not 
contain items with a^a^, &fa% and of erf, a^crf. 

We solved this problem with the finite energy condition §5§ for each sub- 
spaces separately We obtained the conditions for optimal quantum evolution 
in each subspace and calculated the shortest possible time for evolution from 
the initial state to the final one \iPf)- Also, we obtained Hamiltonians 
ffTB]) and ([23D which provide optimal evolution on the subspaces spanned by 
| ft) j I 44) and | f|), | lt)> respectively. 

The Hilbert space of each subspaces is two dimentional and it represented 
by the Bloch sphere. The optimal way of transporting into \ipF) is there- 
fore to rotate the sphere around the axis which pass along the eigenvectors 
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that correspond the largest 00 and the smallest —u eigenvalues of the optimal 
Hamiltonian. Thus, the shortest distance between and is a large 
circle arc on the sphere (geodesic path). 

We used our result for important examples of Uent, Uw^, Uswap, 
U SWAP , UiswAP and U iSWAP gates. In [18] the time of generation of the 
entangler gate is not the shortest possible time. We showed that an en- 
tangler gate can be constructed during the shortest possible time because 
in addition to the Hamiltonian which is considered in [18] our Hamiltonian 
contain J xy and J yx couplings. This allows to reach maximal entangled Bell 
states from the nonentangled ones during the shortest time. Also, we ob- 
tained the conditions for generation of £%]vt' U^swap and U iSWAP gates 
during the shortest possible time. Time of generation Uswap and U qur A P 
t = # is the same as the tune of generation U SWAP in M . This time Ts 
not the shortest possible time t = This is because the generation speed 
v = ^7= is less than maximal value of speed for this Hamiltonian v max = 2u. 
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